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Abstract 

Grover's quantum search algorithm is analyzed for the case in which the initial state is an 
arbitrary pure quantum state \4>) of n qubits. It is shown that the optimal time to perform the 
measurement is independent of |</>), namely, it is identical to the optimal time in the original 
algorithm in which \(f>) = |0), with the same number of marked states, r. The probability of success 
P s is obtained, in terms of the amplitudes of the state \4>), and is shown to be independent of r. A 
class of states, which includes fixed points and cycles of the Grover iteration operator is identified. 
The relevance of these results in the context of using the success probability as an entanglement 
measure is discussed. In particular, the Groverian entanglement measure, previously limited to a 
single marked state, is generalized to the case of several marked states. 
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I. INTRODUCTION 

Grover's quantum search algorithm exemplifies the potential speed-up offered by 

quantum computers. It also provides a laboratory for the analysis of quantum algorithms 
and their implementation. The problem addressed by Grover's algorithm can be viewed 
as trying to find a marked element in an unsorted database of size N. While a classical 
computer would need, on average, N/2 database queries (and N queries in the worst case) 
to solve this problem, a quantum computer using Grover's algorithm, would accomplish the 
same task using merely O(vJV) queries. The importance of Grover's result stems from the 
fact that it proves the enhanced power of quantum computers compared to classical ones 
for a whole class of oracle-based problems, for which the bound on the efficiency of classical 
algorithms is known. Moreover, it was shown |3j that Grover's algorithm is as efficient as 
theoretically possible A variety of applications were developed, in which the algorithm is 
used in the solution of other problems Experimental implementations 

of Grover's algorithm were constructed using nuclear magnetic resonance (NMR) HQ as 
well as an optical device [l^j ]. 

Several generalizations of Grover's original algorithm have been developed. The case in 
which there are several marked states was studied in Ref. It was shown that when 

there are r marked states, Grover's algorithm would be able to find one of them after 
T = 0(yj N/r) queries. A further generalization was obtained by allowing the replacement 



of the Hadamard trans: 
unitary transformation 



orm, used in the original setting, by an arbitrary (but constant) 
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[l7],[lil> as wel] as by the replacement of the n inversion by an 
arbitrary (but constant) phase rotation [l9|. 

Another generalization was obtained by allowing the replacement of the uniform super- 
position of all basis states, used as the initial state of the algorithm in the original setting, 



by an arbitrary pure 
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21 1 or mixed 22 1 quantum state. In this case the probability of 



success of the algorithm, P s , may be reduced, depending on the initial state. Recently it was 
shown that this generalization gives rise to an entanglement monotone that can be used to 
quantify the entanglement in pure states of multiple qubits 23]. 

In this paper we provide a comprehensive analysis of the dynamical behavior and the 
success probability of the quantum search algorithm for arbitrary initial pure quantum 
states. We show that for a given initial state, \(f>), the success probability does not depend 
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on the number of marked states used in the algorithm. We provide an explicit expression 
for P s in terms of the amplitudes of the initial state, and discuss a simple geometrical 
interpretation of it. We then use this approach to calculate P s for states that exhibit various 
symmetries, as well as states that typically appear in quantum algorithms. Two special 
classes of states are identified: fixed points, namely states that remain invariant under the 
Grover iteration as well as two-cycles. The implications of the results in the context of 
entanglement measures for quantum states of multiple qubits are discussed. 

The paper is organized as follows. In Sec. II we briefly describe the algorithm. The 
algorithm in the case of an arbitrary initial quantum state is studied in Sec. Ill, where an 
expression for P s is obtained and shown to be independent of r, and fixed points and cycles 
are identified. The results are discussed in Sec. IV and summarized in Sec. V. 



II. THE QUANTUM SEARCH ALGORITHM 

Consider a search space D containing N elements. We assume, for convenience, that 
N = 2 n , where n is an integer. In this way, we may represent the elements of D using an 
n-qubit register containing the indices, i = 0, ... ,N — 1. We assume that a subset of r 
elements in the search space are marked, that is, they are solutions of the search problem. 
The distinction between the marked and unmarked elements can be expressed by a suitable 
function, / : D — > {0, 1}, such that / = 1 for the marked elements, and / = for the rest. 
The search for a marked element now becomes a search for an element for which / = 1. To 
solve this problem on a classical computer one needs to evaluate / for each element, one 
by one, until a marked state is found. Thus, on average, N/2 evaluations of / are required 
and N in the worst case. For a quantum computer, on which / to be evaluated coherently, 
it was shown that a sequence of unitary operations called Grover's algorithm can locate a 
marked element using only 0(y/N/r) coherent queries of / 

To describe the operation of the quantum search algorithm we first introduce a register, 
\i) = \ii . . .in), of n qubits, and an ancilla qubit, \q), to be used in the computation. We 
also introduce a quantum oracle, a unitary operator O which functions as a black box with 
the ability to recognize solutions to the search problem. The oracle performs the following 
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unitary operation on computational basis states of the register, |z), and the ancilla, \q): 

0\i)\q) = \i)\q®f(i)) (1) 

where © denotes addition modulo 2. The oracle recognizes marked states in the sense that 
if \i) is a marked element of the search space, namely f(i) = 1, the oracle flips the ancilla 
qubit from |0) to |1) and vice versa, while for unmarked states the ancilla is unchanged. The 
ancilla qubit is initially set to the state 

H,= ^(|0>-|1». (2) 
With this choice, the action of the oracle is: 

0\i)\-) = (-l)M\i)\-) q . (3) 

Thus, the only effect of the oracle is to apply a phase of —1 if x is a marked basis state, and 
no phase change if x is unmarked. Since the state of the ancilla does not change, one my 
omit it and write the action of the oracle as 0\x) = (—l)^ x ^\x). 

Grover's search algorithm may be described as follows: Given a black box oracle O, whose 
action is defined by Eq. (JTJ) and a register of n + 1 qubits in the state |0)® n |0)q, the following 
procedure is performed: 

1. Initialization: Apply a Hadamard gate H = ^5(1-1) to each qubit in the register, 
and the gate HX to the ancilla, where X = ( 5 ) * s the not gate, and we write 
matrices with respect to the computational basis (|0), |1)). The resulting state is: 



to>l->* (4) 

where 



N-l 



i=0 



2. Grover Iterations: Repeat the following operation r times (where r is an integer given 
below). 



4 



(a) Apply the oracle, which has the effect of rotating the marked states by a phase of 
7r radians. Since the ancilla is always in the state \—) q the effect of this operation 
may be described by the following unitary operator 

N-l 

^ = E(- 1 ) /W i < x<i. ( 6 ) 

j=0 

acting only on the register. 

(b) (i) apply the Hadamard gate on each qubit in the register; (ii) Rotate the 1 00 . . . 0) 
state of the register by a phase of n radians. This rotation is similar to 2(a), except 
for the fact that here it is performed on a known state. It takes the form 

N-l 

7£ = -2|0><0| + 5»<i|, (7) 

i=0 

where the second term on the right hand side is the identity operator, denoted 
by /; (iii) Apply the Hadamard gate again on each qubit in the register. 

The resulting operation is 



_ H ®njn H m = _ I + 2#® n |O)(O|# 0n = -/ + \v)(v\- 



(8) 



When this operator is applied on the state £V Oj|i) it results in the state 

dj)|i), where a = J2i a i/N- Thus, each amplitude is rotated by n around the 

average of all amplitides of the quantum state. 

The combined operation on the register in one Grover iteration is described by 



U G = -H® n I£H® n I]. 



(9) 



3. Measure the register in the computational basis. 



The optimal number of iterations is 



T = 



-\- 

4V r 



(10) 
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where |_^J is the largest integer which is smaller than x . Moreover, at this optimal 

time a marked state can be found with almost certainty, or more precisely with probability 



* = l-°|-7jj)- (") 



With this performance, Grover's algorithm was found to be optimal in the sense that it is 
as efficient as theoretically possible [4|]. 

Note that the probability P « 1 can be achieved only for specific initial states such 
as the one produced in step 1 of the algorithm above. If this starting state is replaced 



by an arbitrary quantum state, the probability of succes, P s , is reduced 

|2fll2J|. In the 

next Section we analyze the operation of the algorithm in the case that the initial state 
is an arbitrary pure quantum state \<f>). The time evolution of the quantum state during 
the Grover iterations is examined, for a specific choice of the marked states. Special states 
that exhibit fixed points and cycles under Grover iterations are identified. A closed form 
expression for P s is obtained, and a geometrical interpretation for it is presented. It is shown 
that for a given state the success probability does not depend on the number of marked 
states. Finally, some specific quantum states are examined and their utility as initial states 
for Grover's algorithm is obtained. 

III. GROVER'S ALGORITHM WITH AN ARBITRARY PURE INITIAL STATE 
A. The initial state 

Consider an arbitrary pure quantum state \<f>) of n qubits, to be used as the initial state 
in Grover's algorithm. It can be expressed by 

N-l 



5>(0)|i), (12) 



i=0 



where the amplitudes eij(O), i — 0, . . . , N — 1 are complex numbers that satisfy 



N-l 



£>(0)| 2 = 1, (13) 



i=0 



and N = 2 n . The distribution of these amplitudes can be characterized by its moments. 
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The first moment (or average) of the amplitudes is 



— 1 N ^ 

(0) = ^$>(0). 

i=0 



(14) 



The second moment 



N-1 



i=0 



(15) 



satisfies |a(0)| 2 = 1/N, for any state due to the normalization condition (|13j). The 
standard deviation of the amplitude distribution is given by 



Using the equality 



we obtain that 



N-1 



^(0) = ^]Tk(0)-a(0)p 

i=0 



(16) 



<(0) = |a(0)| 2 - a(0) 



(17) 



^(0) 



N 



a(0) 



(18) 



This result can be used to identify two limits. One limit is the equal superposition state, in 
which ai(0) = 1/VN for i = 0, . . . ,N- 1. In this state a(0) = l/v 7 ^ and cr o (0) = 0. In the 



opposite limit there is a large family of states for which a(0) = 0. In this case cr a (0) = l/y/N, 
which is the largest value that cr a (0) can take. 



B. Recursion equations for a pre-defined set of marked states 

Consider Grover's algorithm, searching for one of r marked states where the initial state 
is Denote the set of indices of the marked states by M. The amplitudes of the marked 
states will thus be aj, i G M. In some cases below, we will denote the marked states by 
|mi), . . . , \m r ) and their amplitudes by a mi (0), . . . , a mr (0). The complementary set, which 
includes the indices of the unmarked states, will be denoted by M. 

The time evolution of the amplitudes of the marked and unmarked states during the 
Grover's iterations with an arbitrary pure quantum state as the initial state, was studied in 
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Ref. |2lJ. Starting from the state the amplitudes obtained after t Grover iterations are 
denoted by cij(i), i — 0, . . . , N — 1. 

For a given choice of the marked states, one can consider separately the averages and 
standard deviations of the sets of marked and unmarked states at time t. The averages will 
be 



(t) = ±J>(f) (19) 



r 



for the marked states, and 



°"(*) = w^J2<v ( 2 °) 

for the unmarked states. The standard deviations will be 

ieM 

for the marked states, and 



N — r _ 



for the unmarked states. 

Each Grover iteration consists of two steps. In the first step the phases of all the marked 
amplitudes are totated by n, namely — > — Oj, ieM. In the second step all the amplitudes 
are rotated by 7r around their average, namely dj — > 2a — i = 0, . . . , iV — 1. Using these 
properties, the time dependence was found to be described by the recursion equations 21] 



Oi(t+l) = C(t) +Oi(t) i G M 
ai(t + 1) = C(t) - di(t) i E M 



(23) 



and 



C(t) 



2 

iV 



(N - r)a u (t) - ra m (t) 



(24) 



8 



It was also shown that the standard deviations of the amplitude distributions of the marked 
and unmarked states are constants of motion, namely a m {t) = a m (0) and a u (t) = cr u (0) at 
any time t. 

C. Solution of the recursion equations 

Consider a quantum search with r marked states, using an arbitrary quantum state |</>) 
as the initial state. It was shown in Ref. |2l| that the time evolution of the amplitudes is 
given by 



ai(t) = a m (t) + Adi : jGM 

<n(t) = a u (t) + (-ifAai : i e M, (25) 



where 



N — r 



a m (t) = \ a sin (ut + 5) 



a u (t) = a cos (oot + 5) (26) 



and 



Adi = CLi(0) - a m (0) : i G M 



Adi = Oi(0) - a u (0) : % G M. (27) 
The parameters a and S are given by 



a = Wa„(0) +— a m (0) 

IS — r 



ex P (2^) = ^E^ffi, (28) 
y/N^avW-iy/f a m (0) 

where —tt/2 < Re(<5) < tt/2. Furthermore, it was found that if a measurement is taken after 

t iterations, the probability to find one of the marked states is given by 

P{t) = P -AP cos 2 [ut + Re(8)] (29) 
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where 



p = l-(N- r)a 2 u - l -K 



(30) 



is the highest probability of success that can be achieved for the specific choice of the set 
of marked states denoted by M (if the measurement is taken at the optimal time). The 
parameter K takes the form 



K= (N-r) a u (0) 



_l_ r 



a m (0) - (JV-rK(O) +ra m (0) 



(31) 



The coefficient of the time dependent term is 



AP 



(N-r)a u (0) + T(l m 

(0) 



(32) 



The frequency u given by 



cos uo 



I- 2 -, 

N 



(33) 



or approximately by u = 2^r/N where r <C N. The function P(t) is a sinusoidal function. 
In the analysis of the efficiency of the algorithm we are interested in the largest value Pq that 
this function may reach during its cycle. We are also interested in the number of iterations 
that it would take to reach this maximum for the first time, from the given initial state. 
This optimal number of iterations, for the given set, M of marked states is given by 



1 N /7T 

2V~ V2 



Re{8) 



(34) 



Note that for < Re(5) < tt/2, tm < t namely the optimal number of iterations is equal or 
smaller than for the initial state \rj). On the other hand, for —tt/2 < Re(<5) < the optimal 
number of iterations tm is larger than r. Note that within this analysis, in which the set of 
marked states is assumed to be known, it would be more efficient for —tt/2 < Re(<5) < 0, to 
apply the Grover ierations in reverse order. Of course this does not apply to the actual 
search in which the marked states are unknown. 

There are several issues to consider about the probability Pq. First, the case that we 
considered so far is not the real search. This is due to the fact that in the derivation of 
Eq. (|2H|) it was assumed that the identities of the marked states and their amplitudes are 
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known. In order to obtain the probability of success of an actual search one needs to average 
Eq. over all possible sets of the r marked states. In the case of one marked state this 
yields an average over N possibilities. In general, for r marked states the number of possible 
sets is given by the binomial coefficient ( ^ ) . Averaging over all these possibilities yields 



(a m (0)), which is the average of a m (0) over all possible sets of r marked states. This average 
is used as a small parameter in the series expansion that provides the average of P(t) over 
all these possible choices of the marked states. 

Consider the constant K in Eq. ()31|) . For an initial state \4>) in which all the amplitudes 
are real, there is an exact cancellation and K = 0. Clearly, this simplifies the calculation 



of P significantly. In case that the amplitudes are complex, and the phases of a m (0) and 



a u (0) are different, K may be nonzero. It is easy to see that it is bounded by < K < 



2 



2 



min{2r a m (0) ,2(N — r) a u (0) }. Consider, for example, the case of a single marked state. 
For a specific choice of mi, the constant K can take any value in the range < K < 1. 
However, averaging it over all the N possible choices of the marked state, it is easy to see 
that K is reduced to order 1/N. 

For a given choice of the set of marked states, the optimal time to measure is found to 
depend on the initial state according to Eq. ([34)1 . However, by averaging P(t) [Eq. 
over all possible choices of the set of marked states, Eq. (}3T)j) is replaced by 



P = l-Na 2 a + O[^=), (35) 



and using Eq. (fTHj) we obtain 



P{t) =N a(0) 2 sin 2 (^) + (77=) • ( 36 ) 



Therefore, for any initial state \4>), the optimal number of iterations to be performed before 
taking the measurement remains the same, namely r [Eq. (jlDj)]. The probability of success 
is P s = N 



a(0) 



The averaging process that leads to this result is rather tedius. In the 



next Sections we will present a more elegant approach for the calculation of P s . 
D. Fixed points and cycles 

Typically, the function P(t) is a quasi-periodic function since, generically, the frequency 
a>/(27r) is incommensurate with the sequence of integers counting the Grover iterations. A 
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periodic (rather than quasi-periodic) cycle of the amplitudes is obtained when the frequency 
uj is a rational product of 2n. This condition is satisfied when 



cos 




(37) 



One example of such periodic cycle appears when N/r = 4. In this case uj = n/3, and the 
cycle is of period 6. 

For some initial states and specific choices of the set of marked states, Grover's algorithm 
does not provide any enhancement of the probability to obtain a marked state. Such sit- 
uations appear when AP = 0, thus the success probability P(t) becomes a constant. The 
condition for AP = can be expressed by 



Under this condition the amplitude a = 0, and thus a m (t) = and a u (t) = 0. It is easy 
to see from the geometrical description of the algorithm that in this case the operator Uq 
exhibits cycles of period 2, namely 



The condition for a two-cycle is simply a m (0) = and a u (0) = 0. This is due to the fact that 
under this condition the amplitudes of the marked states are invariant under Uq, while those 
of the unmarked states change sign at each iteration (what matters is the relative phase of 
7i between the marked and unmarked states, while the global phase can be ignored). 
We will also consider fixed points of Uq, namely states that satisfy 



The set of fixed points of Uq consists of two classes: (a) states for which a m (0) = and 
Oj(0) = for each of the unmarked states; (b) states for which a mi (0) = for each of the 
marked states and a u (0) = 0. In the case of a single marked state, \mi), a fixed point is 
obtained when a mi (0) = and a u (0) = (where at least two of the unmarked amplitudes 
are non-zero). 




(38) 



U 2 G \<P) = 10). 



(39) 



U G \4>) = 



(40) 
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Consider the case of two marked states, \ra%) and \rri2). The state 



10) 



1 



(|mi) - \m 2 )) 



(41) 



V2 



is a fixed point of Uq because the amplitudes of all the unmarked states vanish and the 
average amplitude of the marked states vanishes as well. In the case of r marked states, 
I mi), . . . , \m r ), the states 



where a m (0) = 0, are fixed points of Uq- 

We conclude that with the exception of those periodic cycles with periods related to 
the frequency u, the only possible cycles of XJq are the fixed points and cycles of period 2, 
discussed above. 

E. The probability of success for one marked state 

The initialization step of the original algorithm leads to the state |^), given by Eq. 
Consider the case that there are r marked states, |j7ij), i = l,...,r. The operation of 
the algorithm, using the optimal initial state, \r]), and performing the optimal number of 
iterations, r can be expressed by 



We will now examine the case in which the initial state is an arbitrary pure quantum state 
\4>) of n qubits. The purpose of the algorithm is to find one of the marked states, which 
are not known in advance. As we have seen in the previous Section, the optimal number of 
iterations is r, which is independent of the initial state \<p). The final state, just before the 
measurement of the register, will thus be Uq \<f>). We will now calculate the probability that 
a measurement of this state will yield one of the marked states. For simplicity, we will first 
consider the case of one marked state \mx). In this case, for the initial state \rj) we obtain 



r 




(42) 



i=l 




(43) 




(44) 
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For an arbitrary initial state the resulting state will be Uq \<f>), and the probability that 
a measurement of this state will yield the marked state |mi) is 

P s = \{ mi \U^)\ 2 + o(^=^. (45) 
Note that Eq. (J44|) can be written in the form 

<m a |l7S=fa| + oQ=V (46) 
Inserting this equation into ()45|) we obtain that 



\(v\<P)\ 2 + 0(^=), (47) 



namely, P s is given by the overlap of <ft with the equal superposition \rj). 

Note that in the analysis above we used an arbitrary marked state \mi). Since the marked 
state is not known, the calculation of P s should be done by averaging Eq. (|43j) over all the 
N possible choices of the marked state. It will thus take the form 



1 / 1 \ 

where in each of the N terms the operator Uq is constructed such that the marked state is 
the corresponding state It turns out that this averaging is not really necessary, because 
all the N terms in this sum are identical [up to corrections of order 0(l/\/~N)]. 

Eq. (f4Tj) can be used in order to express the success probability P s in terms of the 
amplitudes of \<f>). To this end we will express the states \rj) and \4>) in terms of their 
amplitudes, according to Eqs. (0), and ()12|). respectively, and obtain 



a(0) +0 (49) 



P S = N 

Clearly, P s turns out to depend only on the first moment of the distribution of the amplitudes. 
Initial states that exhibit high values of P s are those for which the average amplitude a is 
large (in absolute value). 
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F. The probability of success for two marked states 



Consider the quantum search algorithm with an initial state |0) and two marked states, 
| mi) and \rri2). The probability that the algorithm will yield one of the marked states is 



P s = Km^Mf + |(m 2 |f^|0)| 2 + O ( -= ) . (50) 



1 

In this case, the operation of the algorithm on the initial state \rj) can be described by 



Uh \v) = ^(Imx) + |m 2 )) + O ( J=) . (51) 



or 



-±=(( mi \ + (m 2 \)U G =( V \ + . (52) 

unlike the case of a single marked state, this equation cannot be applied directly for the 
evaluation of P s . The reason is that P s is determined by the overlaps of U G \<p) with each one 
of the marked states separately and not with their superposition. To overcome this difficulty 
we will use the following identity: 

\{mM)\ 2 + IK|^)| 2 = ~ \({ mi \ + (m 2 |)^)| 2 + ~ \({ mi \ - (m 2 |)^)| 2 , (53) 
where \ip) is an arbitrary quantum state of n qubits. Using this identity we can write 



IKI^I^)! 2 + |(m 2 |^|0)| 2 = i |((mi| + (m 2 |)^|0)| 2 + i |(( mi | - (m 2 |)^|0)| 2 . (54) 

In both terms on the right hand size, we now apply the operator Uq to the left. In the first 
term, according to Eq. this gives rise to the state (77 1. In the second term the state 
on the left hand size of Uq is a fixed point of this operator and thus remains unchanged. 
Therefore 

IKI W>l 2 + IM W>l 2 = l»>l 2 + ~ KM - M)0>l 2 • (55) 

Note that the first term on the right hand size does not depend on the choice of the two 
marked states, similarly to what we obtained for a single marked state. Therefore, the 
averaging over all the possible choices of two marked states out of iV basis states is not 
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necessary. However, the second term does depend on the choise of \mi) and \rri2). Thus the 
averaging is required, resulting in 

N-l 

^s = K^I0)l 2 + iV(iV _ 1) £ km-kmi 2 , (56) 

' mi,m,2=0 

where the cases mi = are excluded from the summation. The first term on the right 
hand side is identical to the expression obtained from P s in the case of a single marked state. 
We will now show that P s for two marked states is the same up to small corrections. To 
this end we will put an upper bound on the value of the second term in Eq. (|56|). using the 
inequality 

|((m 1 |-(m 2 |)0)| 2 <2|(m 1 |0)| 2 + 2|(m 2 |0)| 2 . (57) 

The double summation in the second term of Eq. (|56|) gives rise to sums of the form 
YliLo 1 K^I0)| 2 = 1- The contribution of this term is found to be just a correction of order 
1/N. We thus conclude that the probability of success P s of the quantum search algorithm 
with the initial state |0) and two marked states is identical to Eq. (J47|) . namely, it is the 
same probability as for a single marked state. 

The result of this Section can be easily generalized to the case of r marked states. In this 
case the second term on the right hand side of Eq. ()56|) includes a sum over all the pairs of 
marked states. The conclusion that P s is independent of the number of marked states holds 
as long as r <C N. 

G. Calculation of P s for certain pure initial states 

Consider a state \<p) with amplitudes a^. Since the amplitudes a« = | | , i — 0, . . . , N—l 
are complex numbers and satisfy the normalization condition, the value of \a\ increases as the 
amplitudes become more alligned in the complex plane, namely exhibit a narrow distribution 
of phases 9i, as well as a narrow distribution of the |ai|'s. As the distributions become broader 
the success probability decreases. 

For the state (77), in which all the amplitudes are identical, with a = \/\fN, the success 
probability is P 8 = 1. Consider a state \<p) in which the amplitudes of k of the basis states 
are = 1/Vk, and for all the rest N — k basis states at = 0. For such states the average 
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amplitude is a = Vk/N and thus P s = k/N. Clearly, the success probability increases as 
the amplitude is divided more evenly between the basis states. There are several well known 
quantum states that can now be examined. One of them is the generalized GHZ state of n 
qubits 

|0>ghz = -^(|OO...O> + 111...1)). (58) 

In this state only two of the amplitudes are non-zero. Therefore the success probability 
P s = 2/N — > 0, as the number of qubits increases. A similar situation is encountered for the 
W state of n qubits given by 

\4>) w = 4=(|10 . . . 0) + |010 . . . 0) + • • • + |00 . . . 1)). (59) 
v n 

For this state P s = n/N, which also decays to zero as n is increased. 

A large class of states for which P s vanishes up to 0(l/\/~N) includes all the states for 
which a = 0. This class includes, for example, states in which for each amplitude a« of basis 
state i, there is a state j, with a,j = — aj. In general, a random sampling of states in the 
Hilbert space of n qubits tends to yield states with very small a ~ 1/N. This indicates 
that generic quantum states of n qubits are highly inefficient as initial states for Grover's 
algorithm. 



IV. DISCUSSION 



Recently it was shown that the success probability of Grover's algorithm can be used 
in the construction of an entanglement monotone that quantifies the entanglement of pure 
quantum states of multiple qubits 23]. To quantify the entanglement of the state \(f>), one is 
allowed to perform any set of local unitary operations Ui, U2, ■ ■ ■ , U n on the n qubits before 
feeding the register into the Grover apparatus. These operations are designed to maximize 
the success probability P s of the algorit 1m. For the case of a single marked state, the 



resulting success probability will thus be 



N-l 



^Pmax = max — V |(mi|[/£[/i <g> U- 

Ul,...,Un iV ^— ' 



+ 



mi=Q 



1 



(60) 



17 



where the summation takes care of the averaging over all possible choices of the marked state 
\mi) and the maximization is over all possible sets of local unitary operators Ui, i = 1, . . . , n. 
It was shown that P max can be reduced to the form 

max | (s 

\*)\ + > ( 61 ) 

|si,-.*n> VviV / 

where the maximization now runs over all product states, \s\, . . . , s n ) — |si) ® • • • \s n ), of 
the n qubits. The maximal success probability is thus determined by the maximal overlap 
between \<f>) and any product state of n qubits. The Groverian entanglement measure was 
defined as 



G{(j>) = y/l- P max . (62) 

It was shown that G(<f>) is an entanglement monotone 23l l25l|. This means that for product 
states G(4>) = 0, while for any state \<p) of n qubits it is invariant under any local unitary 
operations on single qubits. Moreover, G(<f>) cannot be increased by any local operations 
on the n qubits (where classical communication is allowed between the parties that perform 
these operations). 

The result presented above, that the success probability P s does not depend on the num- 
ber of marked states, provides a generalization of the Groverian entanglement. It removes 
the restriction that G(<f>) will be defined by a search with a single marked state. The use of 
the Groverian entanglement measure to characterize entangled quantum states of multiple 
qubits may provide useful insight about the nature of these states and their role in quan- 
tum algorithms. To this end one needs to develop efficient computational schemes for the 
maximization over all possible sets of local unitary operations Ui, i = 1, . . . ,n [Eq. ([60)1] . 
or, alternatively, over all sets of product states [Eq. (|BT|)]. 



V. SUMMARY 



We have analyzed the dynamics of Grover's quantum search algorithm for the case in 
which the initial state is an arbitrary pure quantum state \<fi) of n qubits. We have shown 
that the optimal time to perform the measurement that concludes the operation of the 
algorithm is independent of the initial state \<p). It is identical to the optimal measurement 
time of the original algorithm, with the same number of marked states, in which \<fj) = |0). 

18 



An expression for the probability of success P s in terms of the amplitudes of the state \(f>) 
is obtained and is shown to be independent of the number of marked states, r. The fixed 
points and cycles of the Grover operator U G are identified. The relevance of the results in 
the context of using the success probability as an entanglement measure is discussed. In 
particular, the Groverian entanglement measure, previously limited to a single marked state, 
is generalized to the case of several marked states. It is shown that as long as r C A^, G(4>) 
is independent of r. 
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